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For precision coherent measurements with ensembles of quantum spins the relevant Figure-of- 
Merit (FOM) is the product of polarized spin density and coherence lifetime, which is generally 
limited by the dynamics of the spin environment. Here, we apply a coherent spectroscopic technique 
to characterize the dynamics of the composite solid-state spin environment of Nitrogen- Vacancy 
(NV) centers in room temperature diamond. For samples of very different NV densities and impurity 
spin concentrations, we show that NV FOM values can be almost an order of magnitude larger than 
previously achieved in other room-temperature solid-state spin systems, and within an order of 
magnitude of the state-of-the-art atomic system. We also identify a new mechanism for suppression 
of electronic spin bath dynamics in the presence of a nuclear spin bath of sufficient concentration. 
This suppression could inform efforts to further increase the FOM for solid-state spin ensemble 
metrology and collective quantum information processing. 



Understanding and controlling the coherence of quan- 
tum spins in solid-state systems is crucial for preci- 
sion metrology [H-Q, quantum information science |4j- 
@| and basic research on quantum many-body dynam- 
ics 0, 0|- Examples of such systems include Nitrogen- 
Vacancy (NV) color centers in diamond foLndll, phospho- 
rous donors in silicon [12, 13] and quantum dots [14]. A 
Figure-of-Merit (FOM) for precision coherent measure- 
ments with quantum spins is the product of the polarized 
spin density (njvy) and the coherence lifetime (T^), with 
the phase-shift sensitivity S4> scaling as P, [l^, [la| 
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Increasing this FOM requires an understanding of the 
sources of decoherence in the system, and their interplay 
with spin density. For example, in the solid state T2 is 
typically limited by interaction with an environment (i.e., 
bath) of paramagnetic spin impurities; whereas nj^v is 
limited by fabrication issues, the associated creation of 
additional spin impurities in the environment, and the 
ability to polarize (state prepare) the quantum spins. 

The paradigm of a central spin coupled to a spin envi- 
ronment has been studied intensively for many years (see 
; and quantum control methods have been 



17, 1 



e.g 

developed to extend the spin coherence lifetime by re- 
ducing the effective interaction with the environment. In 
particular, dynamical decoupling techniques pioneered in 
the field of NMR ^^^\ have recently been applied suc- 
cessfully to extend the effective T2 of single NV-diamond 
electronic spins by more than an order of magnitude 23 - 
[23. 

Here we study experimentally the FOM of ensembles 
of NV color centers in diamond (Fig. [lja,b)) as a system 



for precision spin metrology such as magnetometry. To 
characterize the dynamics of the composite solid-state 
spin bath that limits the NV electronic spin coherence 
time, consisting of both electronic spin (N) and nuclear 
spin (^"^C) impurities, we apply the technique of spec- 



tral decomposition [26j, |27[ . We describe this technique 
and its experimental application to three different dia- 
mond samples with a wide range of NV densities and 
impurity spin concentrations (see Table 1). With the 
use of dynamical decoupling we realize very large values 
for the spin coherence FOM ^ 2 x 10^'*[ms/cm'^] for two 
of these diamond samples, which is almost an order of 
magnitude larger than in other room-temperature solid- 
state systems [28| ; three orders of magnitude greater than 
has been realized in NV-diamond systems for which T2 
has been measured [H, [2^ ; and comparable to high-NV- 
density samples for which only the inhomogeneous coher- 
ence Ufetime (T2*) was measured |30|. In addition, the 
large NV FOM values realized here are only an order of 
magnitude smaller than the state-of-the-art alkali- vapor 
SERF (Spin Exchange Relaxation Free) magnetometer 
[l6| . which suffers from additional practical limitations 
compared to room-temperature solid-state devices. We 
also find unexpectedly long correlation times for the elec- 
tronic spin baths in two diamond samples with natural 
abundance of ^'^C nuclear spin impurities. We identify a 
new mechanism involving an interplay between the elec- 
tronic and nuclear spin baths that can explain the ob- 
served suppression of electronic spin bath dynamics. In 
future work this effect could be used to engineer dia- 
mond samples for even larger NV spin FOM exceeding 
the current state-of-the-art, e.g., for precision magnetom- 
etry and collective quantum information applications. 

The coherence of a two-level quantum system can be 
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FIG. 1. Quantum spin in a solid-state environment (NV- 
center in diamond) and the applied spin-control pulse- 
sequences, (a) Lattice structure of diamond with a Nitrogen- 
Vacancy (NV) color center, (b) Magnetic environment of the 
NV center electronic spin: ^'^C nuclear spin impurities and N 
electronic spin impurities, (c) Energy-level schematic of the 
NV center, (d) Hahn-echo and multi-pulse (CPMG) spin- 
control sequences. 



related to the magnitude of the off-diagonal elements of 
the system's density matrix. Specifically, we deal here 
with an ensemble of NV electronic spins in a finite ex- 
ternal magnetic field, which can be treated as an effec- 
tive two-level ensemble spin system with quantization 
(z) axis aligned with the nitrogen-vacancy axis. When 
the NV spins are placed into a coherent superposition 
of spin eigenstates, e.g., aligned with the x-axis of the 
Bloch sphere, the measurable spin coherence is given by 
C{t) = Tr[p{t)S^]. Due to coupling of the NV spins to 
their magnetic environment, the coherence is lost over 
time, with the general form 
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where the functional x{t) describes the time dependence 
of the decoherence process. In the presence of a modu- 
lation acting on the NV spins (e.g., a resonant RF pulse 
sequence) , described by a filter function in the frequency 
domain Ft{uj) (see below), it can be shown that the de- 
coherence functional is given bv 3ll 3^ 
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where S{uj) is the spectral function describing the cou- 
pling of the system to the environment. Eq. ([3]) holds 
in the approximation of weak coupling of the NV spins 



to the environment, which is appropriate for the (dom- 
inantly) electronic spin baths of the diamond samples 
discussed here [l^ (see below). 

The spectral decomposition technique allows us to de- 
termine S{lj) from straightforward decoherence measure- 
ments of an ensemble of NV spins, as outlined here. If 
we describe the magnitude of the coupling between the 
NV spins and the spin bath by a fluctuating energy term 
e(i), then S{uj) is the Fourier transform of the second- 
order correlator of e{t): 
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It can be seen from Eq. ([3]) that if an appropri- 
ate modulation is applied to the NV spins such that 
Ft{uj) / (uj'^t) = d{uj — ujq), i.e., if a Dirac (5-function (or 
close approximation) is localized at a desired frequency 
ujq, then x(i) — tS{u!o)/Tr. Therefore, by measuring the 
time-dependence of the spin ensemble's coherence C{t) 
when subjected to such a "spectral ^-function" modula- 
tion, we can extract the spin bath's spectral component 
at frequency u/q- 



S{ujo) = -TTln{C{t))/t. 
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This procedure can then be repeated for different values 
of Wo to provide complete spectral decomposition of the 
spin environment. 

A close approximation to the ideal spectral filter func- 
tion Ft{uj) described above can be provided by a variation 
on the well-known CPMG pulse sequence for dynamical 
decoupling of a spin system from its environment (2ll | 
(Fig. Hfd)). The CPMG pulse sequence is an exten- 
sion of the Hahn-echo sequence [l^ (Fig- Hid)), with n 
equally-spaced 7r-pulses applied to the system after ini- 
tially rotating it into the x-axis with a 7r/2-pulse. The re- 
sulting filter function for this sequence Fcpmg{^) covers 
a narrow frequency region (given by 7r/t, where t is total 
length of the sequence) which is centered at ujq = irn/t, 
and is given by 33[ : 
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We note that the narrow-band feature of the CPMG fil- 
ter essentially defines the effectiveness of this sequence for 
dynamical decoupling. We apply a deconvolution proce- 
dure to correct for deviations of this filter function from 
the ideal Dirac (5-function (see Supplement). 

The composite solid-state spin environment of interest 
here, which is dominated by a bath of fluctuating N elec- 
tronic spin impurities, causes decoherence of the probed 
NV spins through magnetic dipolar interactions. In the 
regime of low external magnetic fields and room tem- 
perature (relevant to the present experiments), the bath 
spins are randomly oriented, and their flip-flops (spin 
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TABLE I. Comparison of key characteristics and extracted 
"average-fit" Lorentzian spin bath parameters for the NV- 
diamond samples studied in this work. 
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Meas. technique 


ensemble 


ensemble 


confocal 


N concentration 


~ 1 ppm 


~ 100 ppm 


~ 50 ppm 


NV density 


~ 10i^[cm-^] 


~ 10^^ [cm-^^] 


~ lO^^fcm"^ 


^•^C concentration 


0.01% 


1.1% 


1.1% 


Echo (1-pulse) T2 


240(6) fis 


2(1) ^s 


5(1) lis 


T2 scaUng 


„0.43(6) 


^0.65(5) 


„0.7(1) 


Max. achieved T2 


2.2 ms 


20Ats 


35/is 


A (expected) 


^ 60 kHz 


^ 6 MHz 


^ 3 MHz 


A (measured) 


30(10) kHz 


7(3) MHz 


1(1) MHz 


Tc (expected) 


~ 15 /is 


^ 0.17 /is 


^ 0.34 lis 


Tc (measured) 


10(5) ^s 


3(2) ^s 


10(5) us 


FOM [ms/cm^] 


2 X 10^* 


2 X 10^* 


10^° 



state exchan ges) can be considered as random uncorre- 
lated events 
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Therefore, the resulting spectrum of 
the bath's coupling to the NV spins can be reasonably 
assumed to be Lorentzian i3ll 13411 : 
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This spin bath spectrum is characterized by two param- 
eters: A is the average coupling strength of the bath 
to the probed NV spins; and Tc is the correlation time 
of the bath spins with each other, which is related to 
their characteristic flip-flop time. In general, the coupling 
strength A is given by the average dipolar interaction en- 
ergy between the bath spins and the NV spins, and the 
correlation time Tc is given by the inverse of the dipo- 
lar interaction energy between neighboring bath-spins. 
Since such spin-spin interactions scale as where r is 
the distance between spins, it is expected that the cou- 
pling strength scales as the bath spin density nspm (i-e., 
A oc Uspin), and the correlation time scales as the inverse 
of this density (i.e., Tc oc \/nspin)- 

Note also that the multi-pulse CPMG sequence used 
in this spectral decomposition technique extends the NV 
spin coherence lifetime by suppressing the time-averaged 
coupling to the fluctuating spin environment. In gen- 
eral, the coherence lifetime T2 increases with the num- 
ber of pulses n employed in the CPMG sequence. For 
a Lorentzian bath, in the limit of very short correlation 
times Tc ^ T2, the sequence is inefficient and T2 oc "nP 
(no improvement with number of pulses) . In the opposite 
limit of very long correlation times Tc ^ T2, the scaling 
is T2 oc n^/^ [H, [22I, (see also recent work on quantum 
dots (si). 

Experimentally, we manipulate the |0)-|1) spin man- 
ifold of the NV triplet electronic ground-state using a 



static magnetic field and microwave pulses; and employ 
a 532 nm laser to initialize and provide optical readout of 
the NV spin states (Fig. ffljc); for details see Supplement 
and 0). 

We applied the spectral decomposition technique de- 
scribed above and extracted the spin bath parameters A 
and Tc as well as the NV spin coherence FOM for three 
diamond samples with differing NV densities and con- 
centrations of electronic and nuclear spin impurities (see 
Table I): (1) An isotopically pure ^^C diamond sample 
grown with chemical vapor deposition (CVD) techniques 
(Element 6), was studied using an NV ensemble micro- 
scope [l5| . This sample has a very low concentration of 
^^C nuclear spin impurities (0.01%), a moderate concen- 
tration of N electronic spin impurities (~ 1 ppm), and a 
moderate NV density (~ 10^'''[cm~^]). (2) A thin-layer 
CVD diamond sample (Apollo) with natural ^'^C con- 
centration (1.1%), high N concentration (~ 100 ppm), 
and large NV density (~ 10^^[cm~'^]) was also studied 
using the NV ensemble microscope. (3) A type lb high- 
pressure high-temperature (HPHT) diamond sample (El- 
ement 6) with natural ^'^C concentration, high N concen- 
tration (~ 50 ppm), and low NV density (~ 10^^[cm~'^]) 
was studied using a confocal apparatus able to measure 
single NV centers. 

In Sample 1 ("'^^C) the overall spin environment for the 
probed NV spins is expected to be simply described by 
a bath of N electronic spin impurities. We studied this 
sample because of its combination of insignificant nuclear 
spin impurities and moderate electronic spin impurity 
concentration, suggesting a favorable (large) FOM value. 
We selected Sample 2 (Apollo) for characterization in or- 
der to compare the NV spin coherence FOM at high N 
concentration and NV density with the FOM at moder- 
ate N concentration and NV density in the ^^C sample. 
Characterization of Sample 3 (HPHT) was important for 
confirming that similar electronic spin bath dynamics are 
determined by the spectral decomposition technique for 
single NV spins and NV ensembles in the high N concen- 
tration limit (even though we expect a small FOM value 
due to the low NV density of the HPHT sample). 

Both the HPHT and Apollo samples have a natural 
1.1% abundance of ^'^C nuclear spin impurities: two or- 
ders of magnitude higher density than in the "'^^C dia- 
mond sample. Nonetheless, for the high N concentration 
in the HPHT and Apollo samples, the simple physical 
picture outlined above, of NV electronic spins coupled 
to independent N electronic and ^"^C nuclear spin baths, 
implies that the spin environment and NV spin decoher- 
encc should still be dominated by the bath of N electron 
spins. 

The results for the ^^C sample are summarized in Fig. 
[21 As displayed in Fig. Efa), we find that the best-fit 
Lorentzian spin bath spectrum (fitted to the average of 
all data points, see Supplement) has coupling strength 
of A = 30 ± 10 kHz and correlation time Tr = 10 ± 5 
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FIG. 2. Application of the spectral decomposition technique 
to an isotopically pure '^^C NV-diamond sample, (a) De- 
rived values for the spin bath spectral functions S„{u}) for 
all CPMG pulse sequences (blue dots) and average values 
at each frequency (red crosses); best fit Lorentzian for the 
mean spectral function {S„{!jj)}n (solid black line); and range 
of best-fit Lorentzians for the individual spectral functions 
5,1(1^) for each CPMG pulse sequence (grey band), (b) Scal- 
ing of T2 with the number n of CPMG pulses: derived from 
NV spin coherence decay data C„{t) (dots); fit of dots to a 
power law T2 = 245/is x n"'*^'®' (solid line); and synthesized 
from the average-fit Lorentzian spin bath spectrum A = 30 
kHz, Tc — 10 /is (open squares). 



/iS, which agrees well with the range of values we find for 
the Lorentzian spin bath spectra Sn{uj) fit to each pulse 
sequence separately (see Supplement), A ~ 30 — 50 kHz, 
Tc ~ 5 — 15 fis. These values are also in reasonable agree- 
ment with the expected "N dominated bath" values for 
A and Tc, given this ^^C NV-diamond sample's estimated 
concentrations of ^^C and N spins (see Table I). 

In Fig. mb) we plot the NV spin coherence lifetime T2, 
determined from the measured coherence decay C„(t) for 
an n pulse CPMG control sequence, as a function of the 
number of CPMG pulses out to n = 128 (dots); and a 
power-law fit to this relationship, from which the scaling 
T2 ex. n'^-^^i^) is found (solid line). This power law ex- 
ponent < 2/3 is consistent with the spin bath Tc being 
about an order of magnitude smaller than the Hahn-echo 
T2 for the ^^C diamond sample (see Table I). The max- 
imal coherence time achieved (using a 512-pulse CPMG 
sequence) was 2.2 ms (see Table I). To check the consis- 
tency of our data analysis procedure and assumptions, 
we used the best-fit Lorentzian spin bath spectrum dis- 



cussed above as an input to calculate the T2{n) scaling 
by synthesizing an NV spin coherence decay C„(<) for 
each CPMG pulse sequence; and then fit to find T2 for 
each synthesized Cn{t) (empty squares in Fig. 2(b)). 
From this process we extract a scaling T2 ^ ■nP-'^^^\ in 
good agreement with the scaling derived directly from 
the Cn{t) data. 



Apollo sample 




10 10 
(0 [rad/s] 

(a) 



^^C Sample: T,o<n° "(^> 




10 20 30 40 50 60 70 80 90 
n (number of CPMG pulses) 

(b) 



FIG. 3. Comparative application of the spectral decomposi- 
tion technique to the ^^C, HPHT and Apollo NV-diamond 
samples, (a) Spin bath spectral functions and associated 
Lorentzian fits using the average-fit method (solid lines) and 
the individual-fit method (color bands), (b) Scaling of T2 
with the number of CPMG pulses: dots - derived from NV 
spin coherence decay data Cn{t); solid lines - fit of dots to 
a power law; open squares - synthesized from the average-fit 
Lorentzian spin bath spectrum. 



The results of the spectral decomposition procedure 
for all three samples are displayed in Fig. [3] and Table ID 
Fig. |3Ua) plots the derived spin bath spectral functions 
Sn(j-^) and associated Lorentzian fits; and Fig. shows 
the measured scaling of T2 with the number of CPMG 
pulses. From the full set of data for each sample we 
extract the coupling strength and correlation time of the 
best average-fit Lorentzian, which we verified to recreate 
correctly the decoherence curves and the T2 scaling, using 
the synthesized data analysis procedure outlined above. 

Importantly, as listed in the bottom row of Table I, 
we found very large values for the spin coherence FOM 
^ 2 X 10^'*[ms/cm^] for the ^^C and Apollo samples. 
These FOM values are 3 orders of magnitude larger than 
previously obtained in other room temperature NV di- 
amond systems for which T2 has been measured, with 
either single spins [291] or ensembles 15[ ; and comparable 
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to the FOM for a high-NV-density sample for which only ing to the expression 38|, 1^ 
T2 was measured [30|. The present results should in- 
form improved precision spin metrology, e.g., magnetom- 
etry using NV diamond. The similar large FOM values 
for the ^^C and Apollo samples are qualitatively consis- 
tent with the trade-off of increased NV density causing 
correspondingly decreased T2 due to increased N impu- 
rity concentration, for a similar N-to-NV conversion effi- 
ciency - 0.1% for these CVD samples. The HPHT sam- 
ple has very low N-to-NV conversion efficiency ^ 10^^, 
and hence has a much smaller FOM. 

Also as seen in Table U there is reasonable agreement 
between the measured and expected values for the cou- 
pling strength A in all three NV-diamond samples, with 
A scaling approximately linearly with the N concentra- 
tion. Similarly, the measured and expected values for the 
correlation time Tc agree well for the ^^C sample with low 
N concentration. However, we find a striking discrep- 
ancy between the measured and expected values of Tc for 
the two samples with 1.1% ^^C concentration and high 
N concentration (HPHT and Apollo): both these sam- 
ples have measured spin bath correlation times that are 
more than an order of magnitude longer than given by 
the above simple electronic spin bath model, though the 
relative values of Tc for the HPHT and Apollo samples 
scale inversely with N concentration, as expected. 

The qualitative agreement between the results of the 
Apollo and HPHT samples indicates that the measured 
behavior of the electronic spin bath at high N concentra- 
tion is not an artifact of the few NVs measured in the 
HPHT sample, and is not related to ensemble averaging 
in the Apollo sample. 

The surprisingly long spin bath correlation times for 
the ^'^C diamond samples, revealed by the spectral de- 
composition technique, imply a suppression of the elec- 
tronic (N) spin bath dynamics at higher nuclear spin 
(^^C) concentration. We explain this suppression as a 
result of the random, relative detuning of electronic spin 
energy levels due to interactions between proximal elec- 
tronic and nuclear spin impurities. The ^^C nuclear spins 
create a local Overhauser field for the N electronic spins, 
which detunes the energy levels of neighboring N spins 
from each other, depending on their (random) relative 
proximity to nearby ^^C spins 37 1. This detuning then 
suppresses the flip-flops of N spins due to energy con- 
servation, effectively increasing the correlation time of 
the electronic spin bath in the presence of a finite con- 
centration of nuclear spin impurities (compared to the 
correlation time in the ^^C sample). 
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The ensemble average effect of the composite 
electronic-nuclear spin bath interaction is to induce an 
inhomogeneous broadening /S.E of the N resonant elec- 
tronic spin transition. The electronic spin flip-flop rate 
R is modifled by this inhomogeneous broadening accord- 



where Ajv is the dipolar interaction between N electronic 
spins. In this physical picture, /S.E is proportional to 
the concentration of ^"^C impurities and to the hyper- 
flne interaction energy between the N electronic and ^"^C 
nuclear spins; whereas A^r is proportional to the N con- 
centration. Taking into account the magnetic moments 
and concentrations of the N electronic and ^'^C nuclear 
spin impurities in our samples, and the large contact hy- 
perflne interaction for N impurities in diamond 4^ 41 1, 
we estimate Ai; ~ 10 MHz and Aat 1 MHz for the 
HPHT and Apollo samples. This estimate implies ap- 
proximately an order of magnitude suppression of R com- 
pared to the bare flip-flop rate {Rtare ~ ^n)j which is 
consistent with the long spin bath correlation times for 
the ^■^C diamond samples, as revealed by the spectral 
decomposition technique (see Table HI. 

In summary, we applied a spectral decomposition 
technique to three NV-diamond samples with different 
composite-spin environments to characterize the bath dy- 
namics and determine the NV spin coherence Figure- 
of-Merit (FOM). Using this technique we realized for 
two samples a FOM ~ 2 x 10"'^^[ms/cm^], which is al- 
most an order of magnitude larger than in other room- 
temperature solid-state spin system; and within one or- 
der of magnitude of the largest spin coherence FOM 
achieved in atomic systems - an alkali vapor SERF mag- 
netometer [l3| (the highest sensitivity magnetometer to 
date). The effective NV spin coherence time achieved 
with the ^^C sample (2.2 ms for 512 CPMG pulses) was 
nearly Ti limited and comparable to the best coherence 
times achieved with single NV spins, suggesting possible 
applications in collective quantum information process- 
ing. For samples with a finite concentration of ^^C nu- 
clear spin impurities, this technique revealed a significant 
suppression of the N electronic spin bath dynamics by in- 
teractions between the nuclear and electronic spin baths. 
Further optimization of the NV spin coherence FOM may 
be possible by engineering this suppression, e.g., with ^'^C 
concentration higher than the natural value. This possi- 
bility is in contrast to the alkali vapor SERF FOM, which 
is fundamentally limited as atomic density is increased by 
spin destruction collisions between atoms [l6| , and there- 
fore holds the promise of NV-diamond devices exceeding 
the current state-of-the-art. 

Finally, the spectral decomposition technique pre- 
sented here, based on well-known pulse sequences and a 
simple reconstruction algorithm, can be applied to other 
composite solid-state spin systems, such as quantum dots 
and phosphorous donors in silicon. Such measurements 
could provide a powerful approach for the study of many- 
body dynamics of complex spin environments, and enable 
identification of optimal parameter regimes for a wide 



6 



range of applications in quantum information science and 

metrology. 
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